Local size segregation in polydisperse hard sphere fluids 
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The structure of polydisperse hard sphere fluids, in the presence of a wall, is studied by the 
Rosenfeld density functional theory. Within this approach, the local excess free energy depends 
on only four combinations of the full set of density flelds. The case of continuous polydispersity 
thereby becomes tractable. We predict, generically, an oscillatory size segregation close to the wall, 
and connect this, by a perturbation theory for narrow distributions, with the reversible work for 
changing the size of one particle in a monodisperse reference fluid. 
PACS numbers: 05.20.Jj,61.20.Gy,64.75.+g 



Understanding the behavior of polydisperse systems is 
relevant to many materials of practical interest. In par- 
ticular, colloidal and/or polymeric fluids generally con- 
tain particles which have, in effect, a continuous distribu- 
tion of sizes (and/or other parameters such as charge and 
chemical composition) . This affects their performance in 
applications ranging from foodstuffs to polymer process- 
ing More fundamentally, colloidal systems also pro- 
vide the closest experimental approach to the 'theorists 
ideal fluid', namely that of perfect hard spheres |^. The 
fact that all colloids are in practice polydisperse (at least 
slightly) must then be taken into account in comparing 
theory with experiment. Only recently has experimental 
work started to clarify in a systematic way the generic 
consequences of polydispersity, such as the partitioning 
of sizes between coexisting phases [|| . 

Despite the continuous interest they have raised, the 
theoretical understanding of polydisperse fluids remains 
far from complete, especially for inhomogeneous cases. 
More is known about partial structure factors in single- 
phase fluids [^ID and liquid-liquid phase equilibrium 
than about crystalline phases g or interfacial properties, 
for example. And, where such inhomogeneous situations 
have been studied Q it has often proved necessary to as- 
sume that only the mean density, and not the size distri- 
bution, can vary in space ||lo[] . This ignores size segrega- 
tion effects, which (globally) influence the phase diagram 
10,^. A similar tendency to local segregation is implicit 
in treatments of binary and ternary hard sphere mixtures 
PI P^ and in polydisperse equilibrium structure factors 
in the homogeneous state 

In what follows, we treat continuous polydispersity 
within a density functional theory (DFT) that properly 
allows for local size segregation. Our work is based on 



a choice of density functional (that of Rosenfeld |14|) 
that has previously been used to study finite mixtures of 
hard spheres. By exploiting the fact that its excess free 
energy density depends on only a small number of lin- 
ear combinations of the particle densities (four 'moment 
densities'), we are able to address the case of continuous 
polydispersity, where the underlying densities are infinite 
in number. This allows us to study, e.g., the effects of 
varying the shape of a smooth size distribution. More- 



over, by a perturbative analysis of the same functional, 
we can distinguish certain generic features that do not 
depend on the shape of the parent, if it remains narrow. 
These two aspects of our work build on recent nonpertur- 
bative |]l5| and perturbative ||] progress in understanding 
polydisperse phase equilibria, extending both approaches 
to the case of inhomogeneous fluids for the flrst time. 

The state of a polydisperse fluid is specified by the local 
number density of each species, p((T, r), with a the parti- 
cle diameter (say). The spatial average of this quantity 
must recover the (known) global size distribution, p(cr), 
which we call the 'parent'. According to density func- 
tional theory, the grand potential is some (unknown) 
functional of p{a,r). Given an approximation to this 
functional, the key problem of polydispersity is the need 
to find, by its minimization, an infinite set of densities 
at each point in space. A naive discretization into (say) 

= 50 species does not make this much easier: the min- 
imization problem remains of very high dimension (we 
have N functions of d spatial coordinates, with d the 
effective space dimension; d ~ \ for a fluid near a flat 
wall) . 

The problem of free energy minimization in a very large 
space arises already in the calculation of phase diagrams, 
where it has been found that, for many purposes, the 
problem can be exactly projected onto a subspace in- 
volving a few linear combinations (or 'moments') of p{a) 
JlSf . The approach requires that the excess free energy 
can be expressed in terms of the chosen moments; this 
is true for many approximate theories including polymer 
mean-field theories and liquid state models such as the 
Percus-Yevick equation of state (PY) (l^Jl^ . 

We now observe that essentially the same simplifi- 
cation is possible with certain (approximate, but well- 
established) density functionals for inhomogeneous liq- 
uids [|l6|. We choose for definiteness that of Rosenfeld 
10 for hard spheres, as subsequently recast by Kierklik 
and Rosinberg fl^ , whose non-ideal part is a functional 
of four 'moment-densities' defined, from the underlying 
density profile p{a, r), as follows (with a = 0, 1, 2, 3): 



TOQ(r) = j dadv' p{a,v')ijJa{(T,\r' — r\ 
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Here (jJq((t, |r' — r|) are four weight functions, selected 
p^ , p^ such that PY is recovered for a homogeneous mix- 
ture. Note that only four moment densities are needed, 
irrespective of the number of components of the mixture; 
this follows from PY itself, whose nonideal part !F'^^{ma) 
involves only the four moments = J da p{a)a°' . How- 
ever, the weight functions Ua are nonlocal: our 'moment 
densities' ma{Y), though intimately related to the four 
PY moments, are not merely local values of them. 
Within this description, the grand potential becomes 

n = j drda {p{a, r) [ln{A'^{c7)p{a, r)) - l] + 

.F- (m„(r)) + {V{a, r) - ^iia))p{a, r)} (2) 

where pier) is the chemical potential of species a in the 
bulk, A(cr) is a thermal wavelength V{a,r) is the 
external potential acting on species a, and T'^^{ma) — 
—mo ln(l — ms) -I- mim2/(l — TO3) -I- m\l (247r(l — m^Y) 
gives the usual PY result [0. Minimization of eq.(||) 
leads to eq.(|^) with 

p(o-, r) = R{a) exp {~[3V{ct, r) + c(cr, r) - c(cr, 00)} (3) 
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Here c((T, 00) is the value of the excess chemical potential, 
and R{(j) the density distribution, in the bulk. 

For a localized potential V (such as a hard wall), and 
for any other case where the densities differ from their 
bulk values only in a finite neighborhood, R{(j) coincides 
with the parent p{a) in the thermodynamic limit [ p^ , 
and is known in advance. Eqs.(^]|^,^) are then closed in 
the low-dimensional function space spanned by {mQ(r)}. 
Their numerical evaluation delivers (via eq.(|2)) the full 
density profile, but nonetheless proceeds as if there were 
only four species present. Note that the four moment- 
density profiles depend (through eqs.(0, ^) on R{a), so 
that the final results, even for the ma{v) themselves, still 
depend on all moments of the parent distribution. 

Of course, not all density functionals in common use 
are of the required 'moment density' form [p^ . But eq.(|^) 
performs as well as most other functionals proposed in 
the literature for example, depletion forces in bi- 
nary fluids are well recovered by it po|| . For numerical 
use below we retain the original weight functions [p^ , 
although some recent modifications are known to give a 
better description of solid phases ; these modifications 
would not alter the conceptual structure of our analysis. 

We now focus on the effect of polydispersity in a fluid 
of hard spheres near a flat hard wall. We fix the parent 
size distribution R{'j) — p{a, 00), and thereby the chem- 
ical potentials of all species in the bulk, and analyse the 
structure as a function of distance z from the wall. We 
have considered fluids with a Gaussian, an exponential. 



and a uniform (top hat) distribution of sizes; each is char- 
acterized by its mean (which we set as ct = 1 without loss 
of generality) and its standard deviation (which is then 
called the polydispersity, p, and is expressed in %). We 
disregard any transverse ordering, which is reasonable for 
a fluid, except possibly close to the fluid/sohd transition. 




FIG. 1. Profiles for the local concentration of several 
species for a flat parent with p = 11.5%, at mean volume 
fraction t; = 0.4. The mean diameter is taken as the length 
unit. The curves for cr = 0.8 and a — 1.2 have been vertically 
displaced to aid comparison. Insets: Size distribution at dif- 
ferent 2. a) Solid line, z = 0.5; dashed line, z = 0.8. b) Solid 
line, z = 1.2; dashed line, z — 1.6. c) Solid line, z = 2.4; 
dashed line, z = 3.2; dotted-dashed line, z = 4.7. 

For all the parents considered, we find significant lo- 
cal size segregation, with a strong, local cross correlation 
between particles of different sizes. This effect can be ob- 
served in fig.0, where we display the local relative concen- 
tration of various species, (l){a, z) = p{a,z)/ J dap{a,z), 
for the case of a top-hat parent with p = 11.5%. Clearly 
visible are strong, ant icor related oscillations in the rel- 
ative amounts of large and small particles, whereas the 
relative concentration of particles close to the mean size 
(cr = 1) is much more nearly constant. Our definition of 
(/)((t) factors out the primary oscillations in the mean den- 
sity (depicted in fig.H below), but we find that the size 
distribution (fig.0, insets) oscillates roughly in quadra- 
ture to the mean density. In this manner, without greatly 
altering the overall density profile, smaller particles are 
accommodated on the inner (near-wall) side of each suc- 
cessive density peak and larger particles on the outer side. 
Thus the density of each species can, within the first few 
layers at least, oscillate with a spatial period close to its 
own diameter. (The same could not continue indefinitely: 
the nature of any size segregation in ordered phases thus 
remains an important, open issue. Moreover, the asymp- 
totic decay of the density profile for any mixture is given 
by a single wavelength l23]). As the polydispersity of the 
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parent increases, the trend we have outhned persists, al- 
though it blurs, and the uniform relative concentration 
of the mean species is lost. The behavior is generic for 
all the different size distributions tried, and can also be 



seen in ternary (or even binary) mixtures |13|. 

We can gain further insight into this behavior by a per- 
turbative analysis for a narrow parent, analogous to that 
used previously for phase coexistence We express 

the particle diameters as cr = a{l + e); we then expand 
in small e around a monodisperse density profile for that 
species. The perturbative result is (to order e) 



' p(™)(oo) ^ ' 

x(l + e[5'(z)-S'(oo)]) 



(5) 



Here p^'"^(z) is the density profile for the monodisperse 
system, and c'(z) is the reversible work coefficient for 
slightly changing the size of one particle, z from the wall, 
within an otherwise monodisperse system. This obeys: 



d 

c'(r) = — c(cr,r)|e=o,p^o 



(6) 



which can be easily computed from eq.(|^). We have as- 
sumed that e ^ 1 for all members of the parent, but not 
that e <IC p; this is why there is no expansion made of 
the factor R{a) in cq.^. Likewise for a hard wall it is 
inappropriate to expand V{<7, z). 



p = 3% 
p = 12% 
p = 19% 
Monodisperse 




FIG. 2. a) Density profiles for the mean species with a flat 
parent, setting (t=1, compared with the profile of a monodis- 
perse fluid at tlie same mean volume fraction, 77 = 0.4. b) 
Reduced local concentration profiles, for a Gaussian parent 
(p = 2.4%) at 77 = 0.45, compared with c' [z) — c'(oo) from 
the monodisperse fiuid at the same -q, thick continuous line. 
The thin continuous line refers to the local concentration of 
the mean species. This curve is not reduced, but has been 
shifted vertically so that it tends to zero at large z. 



These perturbative calculations predict that the profile 
of the 'mean species' is, to order e, identical to a monodis- 
perse fluid at the same overall density. This agrees with 
figj^a, where we compare the mean-species profiles for 
flat parents of varying polydispersity. Curves for other 
parents (not shown), but with matched polydispersity, 
are barely distinguishable on such a plot, even when this 
is over 10%. More generally, eq.(||) shows that all the or- 
der e deviations from monodisperse behavior, for parents 
of different shapes, should fall on a common curve if the 
densities p{a,z) are first normalized by R{a). We have 
verified this by comparing p{a, z)/R{a) for small fixed e, 
among parents of different shapes [|l3| . 

An illuminating application of the perturbation the- 
ory is to the local concentration profiles shown in figj^. 
According to eq.(^, in any region where the external po- 
tential does not depend on the radius of the particles. 



(cr, z) = 4){a, 00) [1 + e (c'(z) - 5'(oo))] 



(7) 



which shows, to order e, that the local concentration of 
the mean species is strictly constant, and that by an ap- 
propriate scaling, the concentration profiles for different 
particle sizes collapse onto a single curve. These predic- 
tions are confirmed from the full (nonperturbative) so- 
lution for a narrow Gaussian parent in fig.|^. The size 
oscillations are, as predicted, directly linked to the re- 
versible work term, c'(z), for changing the size of one 
particle. Although the data collapse is mostly excellent, 
the constancy of the mean species concentration is im- 
perfect close to the wall even for narrow (2.4%) polydis- 
persity. This may be because the perturbative expansion 
itself fails very close to a hard wall, where the poten- 
tial forces the concentration of certain species to be zero. 
(This is visible in the first inset of figj^.) The same prob- 
lem should not arise for a smoothly varying potential. 

We now turn to thermodynamic properties of our sys- 
tem. We have integrated the density and the energy pro- 
files to find the adsorption T and the surface tension . 
In fig.^ we show the F as a function of the volume frac- 
tion, for different parents and polydispersities. It appears 
that, until the parent becomes relatively wide, the values 
of the measured adsorptions do not differ significantly 
from the monodisperse case. Indeed, our perturbative 
calculations show the deviation, and that of the surface 
tension to be of order p^. But in fact the devia- 
tions are only small if, as shown, one uses well-chosen 
moments to scale the plot (effectively, Fto2/too vs 7713). 
These scalings are suggested by (e.g.) scaled particle the- 
ories for the adsorption with different choices, there 
are deviations of up to 30% in the same data. 

In this Letter we have considered the effect of poly- 
dispersity on the structure of an inhomogencous hard 
sphere fiuid. We have shown that the use of moment 
densities makes DFT a practical tool to study inhomoge- 
ncous fluids of continuous polydispersity, in which local 
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size segregation can play a major role. This was illus- 
trated by a study of a polydisperse hard sphere fluid near 
a wall, which clearly shows such effects; spheres segregate 
so as to have density maxima with a period close to their 
own diameter. The generic character of this segregation 
was confirmed by a perturbative analysis. Some of our 
order e size segregation results have analogues in bulk 
phase fractionation of polydisperse fluids. We showed 
that the density of the mean species is unperturbed by 
weak polydispersity; we can also obtain a general rela- 
tion that connects the ratios of differences in the local 
moments (evaluated from p{a, r) at two different spatial 
points) to ratios of higher moments evaluated for the par- 
ent |13|. Such findings generalize similar 0{e) results for 
phase partitioning for these purposes, it is as if, in 
DFT, each point in space counts as a different 'phase'. 
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FIG. 3. Adsorption for different parents. Gaussian: 2.4%, 
7%, 12%; exponential: 9%, 23% (cutoff at a minimum a of 
0.91 and 0.77 respectively), flat: 4%, 12%, 19%. 

We note finally that some of our perturbative results 
(including those just described) are surprisingly general. 
For example, eq.(||), which connects the species densi- 
ties to the local reversible work of enlarging a particle 
in the monodisperse limit, makes no assumption about 
the choice of density functional. Hence it holds for the 
true functional, whatever it may be . Thus the choice 
of physical system (hard spheres), and of approximate 
functional (eq.||) enters only through the particular form 
of c'(r) (eq]|). Eq.(||) is thus an exact result for any 
slightly polydisperse isotropic fluid; moreover the poly- 
disperse feature (a) need not even be size, but could be 
charge, or any other scalar quantity. 
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